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摘要 :【 目 的】 鉴定 中 华 按 蚊 Anopheles sinensis CYP6Y 亚 家 族 基因 ,分 析 它 们 的 结构 和 特征 ,推测 其 可 能 的 功能 。 
【方法 】 以 冈比亚 按 蚊 An. gambiae CYP6Y1 作为 询问 序列 ,通过 双向 Blast 方法 检索 中 华 按 蚊 转 录 组 中 CYP6Y 亚 家 
ee 并 通过 生物 信息 学 方法 分 析 基 因 结 构 特征 及 可 能 的 功能 。[ 结果 】 从 中 华 按 蚊 转 录 组 测序 数据 中 鉴定 出 
2 条 CYP6Y 亚 家 族 基 因 , 分别 命名 为 AsCYP6YI ( GenBank 登录 号 :KF709397) 和 AsCYPGY2 ( GenBank. 登录 号 : 
KF709398) 。 序 列 分 析 显 示 ,AsCYP6Y1 和 AsCYP6Y2 全 长 分 别 为 1713 bp 和 1 815 bp, 分 别 编码 502 和 526 个 氨基 
酸 。 基 因 结 构 分 析 显 示 ,该 亚 家 族 基 因 仅 含有 1 个 相位 1 型 内 含 子 并 与 其 他 P450 基因 形成 保守 的 共 线 性 分 布 。 和 蛋 
白 结构 分 析 显 示 , 这 2 个 基因 编码 的 蛋白 含 P450 特有 的 5 个 特征 序列 和 6 个 底 物 结合 位 点 , 且 均 不 存在 信号 肽 ,其 
亚 细 胞 定位 为 细胞 质 。3D 结构 分 析 显 示 ,AsCYP6Y1 有 18 条 a 螺旋 和 13 股 反 向 平行 的 B dr, AsCYPOY2 有 19 
条 o 螺旋 和 11 股 反 向 平行 的 B 折 苔 。 通 过 同样 的 方法 ,在 达 林 按 蚊 An. darlingi 中 也 鉴定 出 2 个 CYP6Y 亚 家 族 基 
。 系 统 进化 分 析 显示 ,AsCYP6Y1 和 AsCYP6Y2 分 别 与 其 他 3 种 按 蚊 的 CYP6Y1 和 CYP6Y2 聚 成 一 支 , Bootstrap 
值 均 大 于 90% 。 和 替换 率 分 析 显 示 ,中 华 按 蚊 AsCYP6Y1 和 AsCYP6Y2 与 其 他 3 种 按 蚊 同 源 基因 的 Ka/Ks 均 小 于 1。 
相对 进化 速率 分 析 显示 ,中华 按 蚊 CYP6Y 和 CYP6M 亚 家 族 的 相对 进化 速率 均 显 著 快 于 CYP6P 亚 家 族 ,而 CYP6Y 
和 CYP6M 亚 家 族 之 间 没 有 显著 差异 。【 结论 ] 在 中 华 按 蚊 和 达 林 按 蚊 中 存在 2 个 CYP6Y 亚 家 族 基 因 ,之 前 在 冈 比 
亚 按 蚊 和 不 吉 按 蚊 An. funestus 中 也 发 现 2 个 CYP6Y 亚 家 族 基因 ,表明 CYP6Y 亚 家 族 基因 可 能 在 按 蚁 属 广 泛 存 
在 , 且 可 能 为 按 蚊 属 昆虫 所 特有 。 
关键 词 : 中 华 按 蚊 ; CYP6Y 亚 家 族 ; 基因 序列 ; 生物 信息 学 分 析 ; 基因 结构 ; 相对 进化 速率 
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Identification and bioinformatics analysis of genes of the CYP6Y 


subfamily in Anopheles sinensis ( Diptera: Culicidae) 
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Abstract: [Aim] To identify genes of the CYP6Y subfamily in Anopheles sinensis , analyze their structure 
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and characteristics, and deduce their possible functions. [Method] The CYP6Ys cDNA transcripts were 
retrieved and identified by two ways of Blasts using An. gambiae CYP6Y1 as query sequence against the 
transcriptome data of An. sinensis. The structure, characteristics and possible functions of the CYP6Y 
genes identified were analyzed using bioinformatics methods. [ Results] Two genes of the CYP6Y 
subfamily were identified from An. sinensis transcriptome sequencing data, and named as AsCYP6Y1 
( GenBank accession number; KF709397) and AsCYP6Y2 (GenBank accession number; KF709398) , 
respectively. Sequence analysis showed that AsCYP6Y1 and AsCYP6Y2 are 1 713 bp and 1 815 bp in 
length, encoding 502 and 526 amino acids, respectively. Gene structure analysis showed that the genes 
of this subfamily only contain one phase "1" intron and form a conserved synteny with other P450 genes. 
Protein structure prediction showed that the encoded proteins of these two genes contain five P450 
characteristic sequences and six specific substrate binding sites, but have no signal peptide sequence, 
and are localized in cytoplasm. The 3D structural prediction showed that AsCYP6YI has 18 a-helix and 
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13 anti-parallel B-strands, while AsCYP6Y2 has 19 a-helix and 11 anti-parallel g-strands. Two CYP6Y 
genes were also identified in An. darlingi using the same method. Phylogenetic analysis showed that 
AsCYP6Yl1 and AsCYP6Y2 are grouped with CYP6Y1 and CYP6Y2 of other Anopheles species, 
respectively, with the bootstrap values greater than 9096. The Ka/Ks ratio analysis showed that the Ka/ 
Ks ratio of AsCYP6Y1 and AsCYP6Y2 compared with CYP6Y1 and CYP6Y2 of other Anopheles species 
were all less than 1. The relative evolutionary rate test showed that the rates of subfamilies CYP6Y and 
CYP6M were significantly higher than those of the subfamily CYP6P in An. sinensis, while no significant 
difference between subfamilies CYP6Y and CYP6M was observed. [ Conclusion] Two genes of the 
CYP6Y subfamily have been identified and characterized in An. sinensis and An. darlingi, which had 
earlier been found in two other Anopheles species (An. gambiae and An. funestus) , suggesting the wide 
and specific existence of the genes of the CYP6Y subfamily in Anopheles. 

Key words: Anopheles sinensis; CYP6Y subfamily; gene sequence; bioinformatics analysis; gene 


structure; relative evolutionary rate 
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蚊 媒 传染 病 严 重 威 胁 公众 健 康 , Hc PIE f 58 
最 为 严重 ( Rivero et al., 2010; Nkya et al., 2013; 
David et al., 2013) 。 仅 2012 年 ,全球 就 有 2.07 亿 
人 感染 症 疾 ,造成 约 62.7 万 人 死亡 (WHO, 2013), 
喷 酒 杀 虫 剂 以 及 蔗 挂 经 杀 虫 剂 浸泡 的 蚊帐 为 防 控 症 
疾 最 常用 手段 (Enayati and Hemingway, 2010)。 自 
20 世纪 50 年 代 开 始 , 各 类 杀 虫 剂 被 用 于 媒介 蚊虫 
的 控制 (Denholm et al., 2002) ,然而 随 着 杀 虫 剂 的 
广泛 使 用 ,在 全 球 多 达 64 个 国家 已 发 现 蚊 虫 对 杀 贝 
剂 产生 抗 药性 (WHO, 2013 ) 。 鉴 定 蚊虫 抗 药性 相 
关 基 因 ,研究 其 功能 ,阐释 蚊虫 抗 药性 分 子 机 理 对 于 
蚊 媒 疾病 的 有 效 控 制 具有 指导 意义 (Hemingway et 
al., 2002) 。 

蚊虫 抗 药性 机 制 主要 是 靶 标 位 点 突变 (如 
kdr ) 和 解毒 酶 基因 过 量 表达 ( Heckel, 2012) Fé 
虫 体内 解毒 酶 主要 分 为 3 类 :细胞 色素 P450 
(CYP) A BC H BK S- 转 移 酶 和 羧 酸 酯 酶 。 其 中 
P450 为 多 功能 氧化 酶 ,广泛 分 布 于 生物 体 中 ,为 昆 
虫 体内 最 重要 的 解毒 酶 ,并 参与 多 种 生理 活动 ,如 
P450 参与 虹 皮 激素 和 保 幼 激素 的 合成 和 降解 、 植 
物 次 生物 质 代 谢 (Feyereisen，1999 ) 。 与 杀 虫 剂 抗 
性 密切 相关 的 CYP6 家 族 为 昆虫 特有 ,参与 拟 除 虫 
菊 酯 及 有 机 氯 杀 虫 剂 的 代谢 。 经 实验 证 实 冈比亚 
Ti hx Anopheles gambiae CY P6P3 n] 4&6 dg "E cxt o5 RB 
( permethrin ) FIYR 55 5j [li ( deltamethrin ) ; CYP6M2 
能 代谢 3 类 杀 虫 剂 一 一 省 氰 菊 酯 、 氨 基 甲 酸 酯 和 
DDT;CYP6Z1 能 代谢 DDT; 不 吉 按 蚊 An. funestus 
CYP6P9 能 代谢 未 氧 菊 酮 和 澳 氰 菊 酯 ;微小 按 蚊 
An. minimus CYP6AA3 和 CYP6P7 ERIH AE x A28 
Bi URASTE ( cypermethrin ) PI FA BB ( Müller e 
al., 2008; Duangkaew et al., 2011; 
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et al., 2011; Mitchell et al., 2012; David et al., 
2013; Edi et al., 2014) , 

目前 CYP6 Y JV ZA 2E PATE PX EE SE. F AA 
吉 按 蚊 中 有 报道 , 且 仅 发 现 2 个 基因 (Félix et al., 
2010; Irving et al., 2012) ,该 亚 家 族 基 因 与 CYP6M, 
CYP6N, CYP6R, CYP6S 和 CYP6Z 等 亚 家 族 基因 一 
起 构成 一 个 基因 艇 (Nikou et al., 2003; Irving et al., 
2012) 。 仅 在 冈比亚 按 蚊 中 对 该 亚 家 族 基因 的 表达 
做 过 研究 ,通过 分 析 冈 比 亚 按 蚊 感 染 症 原虫 前 后 解 
毒 酶 基因 的 表达 差异 ,发 现 吸食 感染 症 原虫 的 老鼠 
血液 1d 后 CYP6YI f CYP6Y2 分 别 在 脂肪 体 和 中 
肠 中 上 调 表 达 ( Félix et al., 2010)。 

中 华 按 蚊 An. sinensis 是 我 国 症 疾 的 主要 传播 
媒介 ( 刘 起 勇 和 刘 小 波 , 2010; 俞 俊 峰 等 , 2012 ) 。 
本 实验 室 前 期 已 对 不 同 发 育 时 期 中 华 按 蚁 开展 了 转 
录 组 测序 ,本 研究 以 此 为 基础 ,在 转录 组 数据 库 中 搜 
R P450 全 长 序列 ,生物 信息 学 分 析 确 定 其 中 2 条 序 
列 为 CYP6Y 亚 家 族 基因 ,进而 预测 了 已 知 CYP6Y 
亚 家 族 基 因 的 结构 ,特征 和 功能 ,为 进一步 研究 其 功 
能 莫 定 了 基础 。 






































1 材料 与 方法 


1.1 数据 来 源 

中 华 按 蚊 An. sinensis 转录 组 测序 数据 (SRA 登 
录 号 :SRA073189 ) 来 自 于 重庆 师范 大 学 昆虫 与 分 子 
生物 学 研究 所 , P 4e P AE DS H 78 DS V b] EG E f 
An. gambiae 、 不 吉 按 蚊 An. funestus IK PRIR IN An. 
darlingi 和 黑 腹 果 晶 Drosophila melanogaster 等 昆虫 
的 P450 同 源 序 列 下 载 自 NCBI 的 GenBank 数据 库 
(http://www. ncbi. nlm. nih. gov/) 。 
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1.2 序列 鉴定 

以 “P450” 作 为 关键 词 在 中 华 按 蚊 转 录 组 功 
能 注释 结果 中 进行 检索 ,然后 对 得 到 的 预测 序列 
在 NCBI 中 进行 TBLASTX 比 对 。 如 果 与 已 知 P450 
蛋白 序列 一 致 性 高 于 4096 (E-value & e 7) , 则 认为 
其 可 能 是 中 华 按 蚊 的 P450 基因 ( 艾 均 文 ，2008 ) 。 
以 冈比亚 按 蚊 CYP6Y1 ( GenBank. 登录 号 : XPL 
317260.1) 序 列 作为 询问 序列 对 可 能 的 P450 基因 
进行 BLASTX EE xt HX di rp Ae fit CYP6Y 亚 家 族 
基因 。 
1.3 序列 分 析 

使 用 DNAMAN6.0 鉴定 中 华 按 蚊 CYP6Y 亚 家 
族 cDNA 序列 的 开放 阅读 框 并 翻译 成 氨基 酸 序列 ; 
采用 ClustalW 和 GeneDoc 对 中 华 按 蚊 CYP6Y W X 
族 和 其 他 蚊虫 CYP6Y 亚 家 族 进 行 多 重 序列 比 对 
( 姬 继 超 等 ，2013 ) ; 利用 ProtParam 软件 (http:// 
web. expasy. org/protparam/ ) 预测 中 华 按 蚊 CYP6Y 
亚 家 族 蛋 白 的 理论 分 子 量 和 等 电 点 等 ;利用 SMART 
软件 (http://smart. embl-heidelberg. de/) 对 中 华 按 
蚊 CYP6Y 亚 家 族 蛋 白 进行 分 析 ; 利 用 软件 TMHMM 
(http://www. cbs. dtu. dk/services/TMHMM/) 对 中 
华 按 蚊 CYPOY 亚 家 族 蛋 白 进行 跨 膜 区 分 析 ; 使 用 
SignalP 4. 1 $Xf/F ( http://www. cbs. dtu. dk/services/ 
SignalP/) 对 中 华 按 蚊 CYP6Y 亚 家 族 蛋 白 进行 信号 
肽 预测 ; 利用 蛋白 质 亚 细 胞 定位 软件 TargetP 
( http://www. cbs. dtu. dk/services/TargetP/ ) 对 中 华 
fit CYP6Y 亚 家 族 蛋 白 进行 定位 预测 ;利用 NPS 
( http ;//npsa-pbil. ibep. fr/ ) 提供 的 GOR IV 方法 对 
中 华 按 蚊 CYP6Y 亚 家 族 的 蛋白 二 级 结构 进行 在 线 
分 析 ; 先 通过 PSI-BLAST 在 PDB(Protein Data Bank ) 
中 搜索 与 中 华 按 蚊 CYP6Y 家 族 相似 性 高 的 蛋白 再 
采用 SWISS-MODEL ( http://swissmodel. expasy. 
org/ ) 进行 同 源 建 模 及 3D 结构 预测 ,并 使 用 Swiss- 
PdbViewer 4. 1. 0 对 其 构图 。 
1.4 系统 发 育 分 析 

以 中 华 按 蚊 CYP6Y 亚 家 族 序列 作为 询问 序列 ， 
用 BLASTP 软件 在 GenBank 中 搜寻 昆虫 CYP6Y faj 
源 性 序列 并 进行 比 对 ;采用 最 大 似 然 法 以 黑 腹 果 蝇 
CYP6G1 序列 为 外 群 ,用 MEGAS. 1 软件 分 析 中 华 按 
蚊 CYP6Y 亚 家 族 序 列 与 冈比亚 按 蚁 .不 吉 按 蚊 和 达 
林 按 蚊 CYP6Y 亚 家 族 序列 的 系统 发 育 关 系 ,1 000 
次 重复 计算 自 展 分 析 值 。 
1.5 基因 结构 分 析 

利用 在 线 分 析 软 件 FGENESH ( http ;Z/linux1. 





































































































softberry. com/berry. phtml ) 和 GeneWise ( http:// 
www. ebi. ac. uk/Tools/psa/genewise/) 对 NCBI 上 搜 
索 得 到 的 4 种 按 蚊 的 CYP6Y 亚 家 族 基 因 的 基因 组 
DNA 序列 进行 基因 结构 预测 ( Amenya et al., 2005; 
Raghavendra et al., 2012) , 分 析 与 该 亚 家 族 成 复 分 
布 基因 的 基因 结构 ,并 用 Adboe Illustrator CS5 绘制 
其 示意 图 ( Nikou et al., 2003; Irving et al., 2012) 。 
1.6 替换 率 分 析 

先 将 中 华 按 蚊 该 亚 家 族 的 2 个 基因 分 别 与 其 他 
3 种 按 蚊 的 同 源 基因 的 翻译 区 (去 除 终止 密码 子 ) 通 
过 ClustalW 比 对 ,其 结果 用 DnaSP 5. 10. 01 软件 计 
FERRIER XERE (Ka), CERES (Ks) 以 
及 它们 的 比值 (Ka/Ks ) 。 
1.7 相对 进化 速率 分 析 

先 通过 MEGAS. 1 软件 对 中 华 按 蚊 CYP6Y, 
CYP6AA 和 CYP6P 亚 家 族 采 用 邻接 法 构建 系统 进 
化 树 (以 黑 腹 果 晶 CYP6G1 序列 为 外 群 ) ,然后 基于 
其 拓扑 结构 利用 RRTree 软件 计算 相对 进化 速率 指 
标 (mean Ka) 。 






























































2 结果 


2.1 中 华 按 蚊 两 个 CYP6Y 亚 家 族 基 因 的 序列 分 析 
2.1.1 cDNA 序列 分 析 : 通 过 比 对 得 出 中 华 按 蚊 
转录 组 中 ID 号 为 Unigene26541 _5 和 
Unigene27261 5 的 基因 序列 与 冈比亚 按 蚊 
CYP6Yl 基因 的 同 源 性 最 高 ,一 致 性 分 别 为 
84.70% 和 54.50% ,并 且 结 构 域 分 析 显 示 这 两 个 
基因 均 为 P450 超 家 族 基因 ,其 编码 氨基 酸 含有 昆 
虫 P450 蛋白 的 5 个 保守 结构 域 ( 图 1 ~2), 即 
WxxxR, A/GGxE/DTT/S, ExxR, PxxFxPE/DRF, 
FxxGxxxCxG( 血 红 素 结合 区 )( 艾 均 文 等 , 2008 ) 。 
过 分 析 这 两 个 基因 的 cDNA 均 为 全 长 序列 ,将 
Unigene26541 5 命名 为 AsCYP6Yl ( GenBank 登录 
号 :KF709397) ,其 cDNA 全 长 为 1713 bp, 其 中 编 
码 区 1 509 bp ,编码 502 个 氨基 酸 ,其 推导 的 氨基 
酸 序列 如 图 1; 将 Unigene27261 _5 命名 为 
AsCYPG6Y2 ( GenBank 登录 号 :KF709398 ) ,其 cDNA 
全 长 为 1 815 bp, 其 中 编码 区 1.581 bp, 编码 526 
个 氨基 酸 , 其 推导 的 氨基 酸 序列 如 图 2。 

2.1.2  CYP6Yl 和 CYP67Y2 编码 和 蛋白质 的 特性 分 
析 : 初 级 结构 特性 分 析 : ProtParam 软件 分 析 显 示 ， 
AsCYP6Yl 编码 的 蛋白 理论 分 子 量 为 58.2 kD ,等 电 点 
JJ 8. 63 ;AsCYP6Y2 编码 的 蛋白 理论 分 子 量 为 60.4 kD, 
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FxxGxxxCxG 月 








等 电 点 为 8. 74, SignalP 和 TargetP 分 析 显 示 ,2 个 
基因 的 编码 蛋白 均 不 存在 信号 肽 ,其 亚 细 胞 定位 为 
细胞 质 。 和 蛋白 质 结构 域 分 析 表 明 , 中 华 按 蚊 CYP6Y 
亚 家 族 基因 有 5 个 功能 保守 区 :螺旋 C、 螺 旋 1、 螺 旋 
K „Meander 区 和 血红 素 结合 区 (图 3 ) 。 通 过 分 析 发 
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1531 
1621 


图 中 翻译 起 始 密码 子 和 终止 密码 子 被 加 


GCTGTTCGCCGAGCAGTTACAGTTTGGCTTCCAATCCGTTCAGTTCGGAAGC 
[ArgrrGcrTGAGGTAGCAGGAATTGTACTAGCCATCGTGGTGTCGTGTTACGTCTGGCTCAAACGGCGCTACCGGTTTTGGAAGGATCGA 
Mm B WV SACIG IT TWO GE CAS SD Ee uo. Y NOU TZ CK RoCRS Ro SWoAK ADER 
AACGTGGCGTATATCGTGCCCGAGTTTCCGTACGGTAATTTTAAGGCACTCGGTAAGGTGGAACACATCGCACCGATCACGCAGCGGCAC 
Ne Vr A CY ^ Loo BESSER Bn YG SN PK oA C OLoUG kN EB DA CPGE TQ. RUH 
TATGAGCACTTCAAGCGACTGAACGTCCCCTACGGAGGTGTGTTTATGTTGACCTCCCCGCTGCTGTACATCTTCGACACGAAGCTCATC 
TOR UID UESCROCE NOV PAN UG 4 EM DS COT X TL 
AAGACGCTGCTGGTGAAGGACTTCCAGTACTTCCCGAACCGTGGCGTCTACTTCAACGAGCGCGATGATCCACTGTCCGCGCACATGTTC 
K WD TSpSVCKoD CPG WoCBOÓIP. NOR. GS YO TEE NE RD Pb ss MR 
GCGATCGAGGGTCAAAAGTGGCGCGCGCTGCGCGCTAAACTCTCTCCCACGTTCACCTCCGGCCGCATCAAGACGACATTCCCGTTGGTG 
A TB. OAKA WARSA LAR CA OKSCIS op s P TCR AE ese cer TUE qoum SEP. coy 
ACCGAAGTGTGTCGCCGATTCTGTGACCATCTGCTAGAGGAGGTGCAGCGGTCGAACGAGGTGGAGGTGCACGATCTCCTATCGCGCTTC 
THES VOGAR OR F CO OD H: DD Ib UBSE WA SN Er X Boy WW Da RR 
ACCATCGACGTGATCGGCGCGTGCGCCTTCGGCATCGAGTGTAACAGCTTCCGCGAGCCGGACAACGAGTTTCGCCGGTACGGCAAGATC 
T puD Ww LeGTASOGOASGE GO CICER CON CS CER Re ED NE ETR ROY CK 
GCGTTCGATAAGCTGCGCCATTCGCCGCTGGTGGTGTACCTGATGAAGGCGTTCCGGGCGCACGCGAACGCGCTCGGCATGAAGCAGCTG 
ECP D RL OR E S EV OWCYGLI KA RH OA D GOM-K-Q E 
CACGATGACGTGTCGAGCTTCTTCCTGAACGTCGTCAAGGACACGATCGAGTACCGGGAGCGGGAGAAGATCGTGCGCAACGATTTCATG 
Ip DoD esc S qe b WV CK CODO TVL EC YCR E R AERE RN DE HM 





GATCTGCTGGTGAAGCTTAAGAACACGGGGCGGCTGGAGGACGGCGGGGAAGAGATCGGCCGGTTGACCTTCGAAGAAATTGCCGCCCAG 
DAT LONE ERAN GC Eb DG GTE E DIG RS LoT ERESCE Th Arg 
GCGTTCATATTCTTTACCGCCGGCTACGACACCTCGTCAACTGCCATGTCGTACACGCTGTACGAGTTGGCGCTGAACCAGGCGGCGCAG 





ALI ET wh GD TT Sx»): DoOASM 4T I Ys ES4D OA OLOUNCQ^ ACA CQ 
GAAAAGGCTCGCCAGTGTGTCAAGGACACGCTGCAGAAGTACGGCGGGCAGGTTACGTACGAGGCCACTCACGATATGCCCTATCTCGAG 
ECK ACRO (CV D DL A £X o 60-—G.Q VUDuWX«EUERCDIOHBODVMCPBÓOYCL-E 
CAGTGTATCAGTGAAACACTCCGCAAACATCCACCAGTGGCGATTCTGGAGCGCAATGCGGACAA AGACTA TCGCCTACCAGACAGTGGG 
QC. TIS "E T-CL-ROK*H CP BOX SGASCD DOES RON MA DK DY LR mE DS 6 
TTGCTGCTACGCCGTGGCCAGAAGATCATGGTCCCGATCTACGCGATGCACCGCGATCCGGAACACTTTCCCGAACCGGACGAGTACCGG 
L-—LLE RARO Gf0.K cE VB MR DP mb HOPES DB: YR 
CCGGAACGGTTCGATCCGGAAGAGGTAGCTAAGCGGGACTCGCACTGCTACCTTCCGTTCGGTGAGGGACCACGCATATGTATCGGCATG 
POBECRCBES De Er V vAr Be LP Re Ts OP DoG-M 
CGGTTCGGATCCATACAAGCGAAGGTAGGCCTGGCCAATCTGCTGAATCGCTTCCGGTTCAGTGTGTGTGATCGGACGCAGATTCCGGTG 





ES 
CAGTACTCGAGGACGAACTTCATTCTCGGCCCGGCCAATGGTGTTTGGCTTOGTGCTGAGACTCTQTAQGGTCTCCAATAAGTTAAGAAG 
OCT EOROTION E b GP MA GO OE CLGOESA ROT OL o* 
TAGTCAATTGCAACACGTTTGCATGTTAGTTTCCAAATGGTGCTTTACTGAATAAAAGAATGTAAATTTTATTCTCATCCAAAGCATGCT 
GTGGGGGGACGTGTGTCGTCCATTTCTTTACAAATAAACGA 
图 1 中 华 按 蚊 AsCYPOYI cDNA 及 其 推导 的 氨基 酸 序列 
Fig.1 cDNA and deduced amino acid sequence of AsCYP6YI from Anopheles sinensis 











E AE BEC CAE RR Ah E 751 , P450 保守 结构 域 WxxxR, A/GGxE/DTT/S, ExxR, PxxFxPE/DRF 和 
日 粗 体 加 下 划 线 标 出 。The start and stop codon are boxed, the numbers on the left are the positions of nucleotides and amino acids on the 
sequences, and the insect P450 conserved motifs including WxxxR, A/GGxE/DTT/S, ExxR, PxxFxPE/DRF and FxxGxxxCxG are underlined in bold. 
图 2 [i] The same for Fig. 2. 


Jil AsCYP6Y1 和 AsCYP6Y2 的 5 个 保守 结构 与 昆虫 
P450 的 保守 结构 还 是 存在 一 定 差异 : CYP6Y 
Meander [X ( PxxYxPE/DRF ) 中 的 第 4 个 氨基 酸 由 栈 
氨 酸 代替 了 昆虫 P450 ( PxxFxPE/DRF) F WEN A 
酸 ,两 者 仅 在 葵 环 上 有 1 个 羟基 的 区 别 , 均 为 芳香 族 
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氨基 酸 。6 个 底 物 结合 位 点 在 图 3 中 用 红 框 标 出 ,6 


个 位 点 


点 的 保守 性 都 很 高 。 


CTGAGTGCCCTTTCCGCT 
ATCGCCACCGCAGATAGCAATTGTGTCGTGTGTGTGTGTACTCCGCGTTCTTCTACAACCACGGTGTGGCCCGTTAGTAGTCGCTCGGCG 
[ATdGAGCTCGCGGTAGCAGTGGCCAC CTATCTGGTGGCGTTCGTAGCGACGCTCTGTCTGGCCGTGTATCTGTACGTCGGGCGTAAGTAC 

机 
TCCTTCTGGAAGCGGCACCATGTACCGTACGTCCAGCCGGAACTACCGTTCGGGAACTTTAAGGAGATGGGCAAATCGATCCACCCTGCG 
SABE VERR HCHP PEEN X QuPCEOLSP DG NCESUECEIMGUK SA THEE A 
CATCTCAGCAAGCGGTTCTACGACGAGTTCAAGGGAAAGGGCCGTGGAGGCTTCGTTGGGTTGTACATCTTCCTGCAGCCAGTGTTGTTG 
HD o5 RE E GEE AD EF RO E A Gp Ge CY S a e S a a A T 
GTGACCGATTTGCGGCTAGCGAAGCGGATCCTTATCGAAGACTTTCACCACTTCCCCGACCGAGGCGTTTACTTCAACGAGAAGGACGAT 
了 
CCACTGTCGGCGCATCTGTTCTCGCTCGAGGGCAAACGGTGGAAGGATCTGCGGGGGAAAATCAGTCCCACGTTCACCAGCGGGCGTATG 
PARE ST A SH CE B SAL CEXGCROR WR DUE RERE 0E Vo cPAED B STi scG RM 
AAGGCGGCCTTTCCGATGGTATTGAAAATCGGTGAACAGTTCTGCGAGTACCTGCGGGACGAGTGCTCCGGTGGGAAGGTGTTGGAAGTT 
K-A Bb PP MN T KS Deco EL DA E76. Sneek SEES 
CGCGATCTCATGGCACGGTTTACGACCGACGTAATCGGGACGTACGCGTTCGGGTTGGAGTGCAACAGTTTCCGGGATCCGCAGAACGAA 
Ro Die AM AUGE T ED SIOC-T RG bs QS UN GE. 
TTCCGGCGCATCGGACGCAAACACTTCGATACGCCGCGCAACCATCCGCTGAAGGTGTTCGTGATGAAGACGTTCCGTGGGTTGGCGAAC 
FOROR G CIRKCK- HA PADOT «POOR AN CHGB EARE W E VA R T ER G La AN 





GCGATGGGGCTGAAGCTGCTGCACGACGATGTGGCCGCATTCTTTAGCGGTGTCATCCGGCAGACTGTCGCCTACCGGGAGCGACACGAA 
A M-5 LOK L LH. D.D V A.A b P.5 6 V I EQ VA RE RH 
ACCGAGCGGCGGAACGATTTTCTTGACCTGCTGATAAGCCTGAAGAACACGGGCATGCTGGAGAGCTCGGGCGAGATAATAGGCCGCCTG 








TAIE iR RE SNC CODE HC E DS ECESSLSS * Bu K NS UE GE 5G HBS P OD RO 
ACGCTGGACGAAATTGCCGCCCAAGCGTTCATCTTCTTTACGGCCGGCTTCGAGACGTCGTCGAGTGCGATGACCTACACGCTGTACGAG 





Ipse CES ASOA CQ ATE SE CTEUBC CE A GP -ES S DAUM STEEL TE 
ATGGCACTGAACGAGGAGGTCCAGCGGAAGGCACGGGACTGCGTGATGGATGCGCTCGAGAAGCATGACGGCAGCATGACGTACGAGTCG 
MA N ESES Vn a RARE A RODOC V OMODeACL. EOK HD GS MET Y-E $ 
TTCAAAAACATGCTCTATCTGGACCAGTGCATCTATGAAACGCTCCGCAAATATCCACCGGTGGCCATCCTCGAGCGAACCGTCGCGAAG 
Esp NAMA SY TES OQ SUE SY IBSSRSECRIRCYUBEPBSVUASMLIE BOROTOAWAUKR 
CCCTACAAAATCCCCGACACAAACGTTACCCTGCATCGTGGTATGAAGATCATGATCCCCGCGTACGCCATCCACCACGACGAGGAGTTC 
Posee Sp2 Do STA UNS CVO OEC TS CHOIR GS MAR MALS Bo ho N v S EUCH OH DE EE 
TACCCGGACCCGTCGACGTACGATCCGGGTCGCTTTACCCCAGACCGGATGGCCCGGAGGGACCCGTGCGCCTACCTACCGTTCGGCGAA 


yb DP ST To DE GC RE T P Do RM Au RR D PO ÀA^Y Ps Pe 





GGTCCCCGGATATGCATCGGCCTGCGGTTCGGTATGATGCAGGCGCGGATCGGTTTGGCATTACTGCTGAAGCACTTCAAGGTGGTCCAG 
Gd ALSCECDCK mem CA v 
TGTGATGAGACGGACGTGCCCCTCACCTATTCCACGCGAGCATTCGTACTGACGCCCATCAATGGGGTTCGATTGAGGTTGGAGAAATAC 


GSP-CR-I-C L-G LRE GIM M Q A.R 





CU DIES T: Doc BAD ro RS AE OY e ET. SEAN GN. (VEG IZ CR TE E: E N 
GAGCACGTGCCTGGAATTGGAATGGCCTCCGATAAGTTGCTGGATAT TTAAATGCATCAATTAGCCACCGAATAATTCOGTTACCTGCAA 
EC HU WP eG DG M Aus DSK: TED CD ox 
TGGACCGGAATAGTTGTCTTTCATAGCGCTCCAACAAATCACAGAAAATTAAATGTTTTTTTTTTTCTTATCTTCGTAATTCTGAAA 


图 2 ”中华 按 蚊 AsCYP6Y2 cDNA 及 其 推导 的 氨基 酸 序 列 
Fig.2 cDNA and deduced amino acid sequence of AsCYP6Y2 from Anopheles sinensis 





33.2796 和 15. 94% ;在 AsCYP6Y2 的 二 级 结构 中 ， 
无 规则 卷曲 占 的 比例 最 高 为 46. 58% , a- 螺旋 和 延 


蛋白 质 的 空间 结构 预测 :应 用 GOR IV. 方法 对 
中 华 按 蚁 CYP6Y 亚 家 族 和 蛋白 二 级 结构 进行 在 线 分 
析 , 在 AsCYP6Y1 的 二 级 结构 中 ,无 规则 卷曲 占 的 比 
例 最 高 , 达 50. 8096 , a- 螺 旋 和 延伸 链 分 别 占 


伸 链 分 别 占 39.35% 和 14.07%。 由 此 可 推测 ,无 
规则 卷曲 是 AsCYP6Y1 和 AsCYP6Y2 二 级 结构 中 最 
大 量 的 结构 元 件 ,a- 螺 旋 和 延伸 链 分 散 于 整个 蛋白 质 
中 。 通 过 PSI-BLAST 的 搜索 ,AsCYP6Y1 和 AsCYP6Y2 
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3 中华 按 蚊 CYP6Y 3 




















其他 按 蚊 CYP6Y 亚 家 族 基因 氨基 酸 序列 比较 


Fig.3 Amino acid sequence alignment of AsCYP6Y gene from Anopheles sinensis with CYP6Y genes from other Anopheles species 
Ad: 达 林 按 蚊 An. darling; Af; 不 吉 按 蚊 An. funestus; Ag: DJ EGET An. gambiae; As; "BET An. sinensis. 图 中 用 线 框 标 出 的 为 P450 的 

















5 个 保守 结构 域 ,用 方 框 标 出 的 为 6 个 底 物 结合 位 点 ,黑色 、 灰 色 和 白色 阴影 分 别 表示 氨基 酸 序列 保守 性 为 100% ,80% 和 8096 以 下 。Five 


conserved domains are line-boxed out in the figure. Six substrate recognition sites ( SRSs ) are marked with rectangles. Black, grey and white shade 


denote amino acids with 10096 , 80% and below 8096 identity, respectively. 


与 人 CYP3A4 氮 基 酸 序列 一 致 性 分 别 为 32. 096 和 
30. 0% , 被 选 作 同 源 建 模 的 模板 。 通 过 SWISS- 
MODEL 建 模 预 测 得 到 AsCYP6Y1 和 AsCYP6Y2 的 
蛋白质 三 级 结构 。 通 过 Swiss-PdbViewer4. 1. 0 对 其 
构图 如 图 4 所 示 , 结果 分 析 显 示 : AsCYPGYI 有 18 
条 a 螺旋 和 13 股 反 癌 平行 的 B HA, AsSCYPOY2 有 
19 条 a 螺旋 和 11 股 反 向 平行 的 B Dr. 
2.2 CYP6Y 亚 家 族 基 因 系 统 发 育 分 析 

氨基 酸 多 重 序 列 比 对 分 析 显 示 , AsCYP6Y1 与 
AfCYP6Y1 一 致 性 最 高 , 达 84. 7% , AsCYP6Y2 与 
AfCYP6Y2 一 致 性 最 高 , 达 78. 2% 。 选 择 中 华 按 蚊 
































CYP6Y1 和 CYP6Y2 及 其 他 按 蚊 同 源 氨 基 酸 序列 ， 
RERE CYP6G1 为 外 群 构建 系统 发 育 树 ,结果 
表明 AsCYP6Y1 和 AsCYP6Y2 分 别 与 其 他 3 种 按 蚁 
的 CYP6Y1 和 CYP6Y2 聚 成 一 支 , bootstrap 值 均 大 
于 90% (图 5)。 
2.3  CYP6Y 亚 家 族 基因 结构 分 析 

基因 结构 分 析 预 测 表明 AsCYP6Y 亚 家 族 基 因 仅 
含有 1 个 内 含 子 ,该 内 含 子 将 保守 结构 螺旋 
(ETLR) 中 编码 谷 氨 酸 的 密码 子 分 隔 开 , 且 都 为 相位 
1 型 内 含 子 ( 内 含 子 位 于 密码 子 内 的 第 1 NFR 
之 后 ) ,这 与 之 前 对 冈比亚 按 蚊 和 不 吉 按 蚊 中 CYPOYI 
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图 4 中 华 按 蚊 CYP6Y1(A) 和 CYP6Y2 
(B) 的 三 维 结构 预测 

Fig.4 Predicted 3D structure of CYP6Y1 (A) and 

CYP6Y2 (B) from Anopheles sinensis 








的 分 析 结 果 相 同 (Amenya et al., 2005) , CYP6Y1 的 
内 含 子 、 外 显 子 1、 外 显 子 2 长 度 范围 为 分 别 为 76 ~ 
95,1 093 ~1 096 和 416 ~422 bp, CYP6Y2 BS 
子 、 外 显 子 1、 外 显 子 2 长 度 范围 为 分 别 为 63 -93, 
1 114 ~1 123 和 419 ~464 bp( 表 1) ,中 华 按 蚊 中 该 
亚 家 族 基 因 与 AsCYP6N2,AsCYP6R1 和 AsCYPGSI 
成 簇 分 布 。 该 亚 家 族 基因 及 与 其 成 簇 分 布 的 几 个 基 
因 内 含 子 分 布 如 图 6, 从 图 中 可 以 看 出 4 种 按 蚊 中 
该 亚 家 族 基因 周围 基因 分 布 相同 ,存在 保守 的 共 线 
性 分 布 。 
2.4  CYP6Y 亚 家 族 基 因 的 替换 率 分 析 

中 华 按 蚊 AsSCYPOYI $1 AsCYPGY2 与 3 种 其 他 按 蚊 
同 源 基因 的 Ka 均 大 于 0.85 ,Ks 均 小 于 0.20 ,其 Ka/Ks 
值 如 表 2 所 示 ,Ka/ Ks 值 均 小 于 1, 表现 出 纯化 ( 负 ) 选 
择 。 用 于 分 析 基 因 的 GenBank 登录 号 见 表 1 。 

95 [— — AfCYP6Y1 AFM08404.1 
91 ———— AgCYP6Y1 XP 317260.1 

















100 L— — AsCYP6Y1 KF709397 


AdCYP6Y1 EFR25763.1 
AdCYP6Y2 EFR25764.1 
AsCYPY2 KF709398 
99 AfCYP6Y2 AFMOS405.1 
100 —— AgCYP6Y2 XP 317261.1 





0.2 


DmCYO6G1 ADJ57460.1 


图 5 4 种 按 蚊 CYP6Y 亚 家 族 基 因 基 于 最 大 似 然 法 的 系统 发 育 关系 

Fig.5 The phylogenetic relationship of CYP6Ys of four Anopheles mosquito species based on Maximum Likelihood method 
各 序列 名 称 后 为 序列 的 GenBank 登录 号 ;使 用 黑 腹 果 蝇 Drosophila melanogaster ( Dm) ff] CYP6G1 序列 作为 外 群 ;大 于 50% 的 bootstrap 值 标记 在 
树 的 分 支 节点 上 ;标尺 代表 系统 发 育 距 离 。The GenBank accession numbers are listed after the sequence names. The homologous CYP6G1 from 


Drosophila melanogaster ( Dm) is used as the outgroup, percentage bootstrap values higher than 50% are indicated on each branch, and the scale bar 


indicates the phylogenetic distance. 











R1 4 种 按 蚊 CYP6Y 亚 家 族 基 因 内 含 子 形式 和 外 显 子 长 度 


Table1 The pattern of introns and the length of exons among CYP6Ys of four Anopheles mosquito species 





基因 名 GenBank 登录 号 内 含 子 长 度 (bp) 外 显 子 1 长 度 (bp) 外 显 子 2 长 度 (pp) 
Gene name GenBank accession no. Intron size Exon 1 size Exon 2. size 
AsCYP6YI ATLV01014697. 1 85 1 093 416 
AgCYP6Y1 AAAB01008964. 1 76 1 093 422 
AfCYP6YI JQ711183. 1 76 1 093 416 
AdCYP6Y1 ADMHO1000000 82 1 096 416 
AsCYPGY2 ATLV01014697. 11 67 1117 464 
AgCYP6Y2 AAAB01008964. 1 65 1117 443 
AfCYP6Y2 JQ711183. 1 63 1114 431 
AdCYP6Y2 ADMHO1000000 93 1.123 419 
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Fig.6 The genomic distribution of CYP6Y subfamily genes from four Anopheles mosquito species 





种 名 之 后 为 CYP6Ys 所 在 scaffold 的 GenBank 登录 号 ;黑色 和 矩形 为 外 显 子 ,无 填 矩 形 为 内 含 子 ,垂直 箭头 表示 基因 间 区 域 ,水 平 箭头 表示 每 个 基 
因 的 5' -3“ 方 向。The GenBank accession numbers of scaffolds that those genes are located on are listed after the species names. The black-filled 






































rectangles denote protein-coding sequences , and unfilled rectangles denote introns, vertical arrows denote intergenic regions while horizontal arrows denote 


the 5' —3' orientation of each gene. 


x2  AsCYP6YI1 和 AsCYP6Y2 与 3 种 其 他 按 蚊 
同 源 序列 的 Ka/Ks 值 
Table 2 The Ka/Ks values of AsSCYP6Y1 and AsCYP6Y2 


with their homologues of other three Anopheles species 





2.5 ”中华 按 蚊 CYP6Y 亚 家 族 基 因 的 相对 进化 速 
率 分 析 

相对 进化 速率 分 析 显 示 , 中华 按 蚊 CYP6Y 亚 家 
族 基 因 的 进化 速率 显著 快 于 CYP6P 亚 家 族 (mean 
Ka, 0. 728 vs. 0.648, P =0. 00968) ,CYP6M 亚 家 族 
基因 也 显著 快 于 CYP6P (mean Ka, 0. 708 vs. 
0.648, P 20.0462) ,而 CYP6Y 亚 家 族 和 CYP6M 亚 
家 族 之 间 没 有 显著 差异 (mean Ka, 0. 728 vs. 
0.708, P=0.524)( 图 7)。 
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AsCYP6M3 ATLVO1012548.1 
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AsCYP6Y2 KF709398 
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AsCYP6P4 ATLV01011656.1 
AsCYP6P3 ATLV01011653.1 
AsCYP6P5 ATLV01011654.1 
DmCYP6G1 ADJ57460.1 








AsCYP6YI| AsCYP6Y2 

AdCYP6Y1 0.112 
AfCYP6YI 0.075 
AgCYP6Y1 0.091 
AdCYP6Y2 0. 181 
AfCYP6Y2 0.105 
AgCYP6Y2 0.162 

0.708 

0.728 
0.648 
图 7 ”中华 按 蚊 CYP6Y 亚 家 族 基因 基于 邻接 法 的 相对 进化 速率 





Fig.7 The relative evolutionary rates of the CYP6Y genes from Anopheles sinensis based on Neighbor-Joining method 
各 序列 名 称 后 为 序列 的 GenBank 登录 号 ;使 用 黑 腹 果 晶 CYP6GI 序列 作为 外 群 ;CYP6Y ,CYP6M 和 CYP6P 亚 家 族 的 Ka 值 标记 在 树 的 分 支 节 点 
Eo The GenBank accession numbers are listed after the sequence names. The homologous CYP6G1 of Drosophila melanogaster is used as the outgroup, 


and the Ka values of CYP6Y , CYP6M and CYP6P subfamily are indicated on each branch. 





3 讨论 





CYP6Y 亚 家 族 基因 仅 在 冈比亚 按 蚁 和 不 吉 按 
蚊 有 报道 ,并 仅 有 2 个 基因 成 员 ( Felix et al., 2010; 








Irving et al., 2012) 。 本 研究 表明 ,中 华 按 蚊 CYP6Y 
亚 家 族 也 只 有 2 个 基因 (图 1 和 2) , 达 林 按 蚊 基 因 
组 中 CYP6Y 亚 家 族 也 只 有 2 个 基因 (图 6)。 本 研 
究 不 仅 增 加 了 CYP6Y 亚 家 族 基 因 序 列 , 同时 表明 
CYP6Y 亚 家 族 基 因 可 能 是 按 蚊 属 昆虫 所 特有 。 
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中 华 按 蚊 AsCYP6Y1 和 AsCYP6Y2 均 含 有 
P450 特有 的 5 个 保守 区 和 6 个 底 物 结合 位 点 。 
P450 亚 家 族 成 员 间 的 一 致 性 大 于 55% (Nelson et 
al., 1996; Feyereisen , 1999) ,但 中 华 按 蚊 CYP6Y 亚 
家 族 的 这 2 个 基因 间 的 一 致 性 为 52. 5%, 低 于 
55% , TEI EPR EY 不 吉 按 蚊 和 达 林 按 蚊 中 该 亚 家 
族 基因 间 的 一 致 性 均 低 于 5596 ,但 高 于 46% ,这 与 
哺乳 动物 CYP17 的 研究 结果 (Nelson et al., 1996) 
相似 。 

基因 翻译 区 的 Ka/Ks 值 可 以 反映 基因 承受 的 
选择 压力 ,如 Ka/Ks <1, 则 认为 有 纯化 选择 的 压力 ， 
Ka/Ks 值 越 小 表明 选择 压力 越 大 , 非 同 义 替 换 可 能 
对 其 产生 有 害 作 用 ,不 利于 其 生存 (Wagner，2002; 
周 琦 和 王 文 ， 2004)。 中 华 按 蚊 ASCYPOYI 和 
AsCYPGY2 5j 3 种 其 他 按 蚊 同 源 基 因 Ka/ Ks 值 均 小 
于 1, 介 于 0.07 ~0.20 之 间 , 表 明 该 亚 家 族 基因 可 
能 纯化 选择 压力 大 ,暗示 其 对 于 中 华 按 蚊 的 生存 起 
着 重要 作用 。 

ERER EREI H P450 基因 是 研究 其 起 
源 和 进化 的 关键 (Feyereisen ,2011 ) , CYP6 家 族 基 
因 多 成 复 分 布 ,在 冈比亚 按 蚊 和 不 吉 按 蚊 中 该 家 族 
基因 均 构 成 了 3 个 基因 艇 (Nikou et al., 2003; Irving 
et al., 2012) ,中 华 按 蚊 基 因 组 数据 分 析 显 示 CYP6Y 
亚 家 族 基 因 与 另外 3 个 CYP6 家 族 基因 构成 了 1 个 
基因 簇 , 图 6 显示 4 种 按 蚊 中 CYP651 与 其 他 4 个 
基因 的 转录 方向 相反 ,并 且 除 达 林 按 蚊 的 CYP6S1 
外 的 其 他 基因 内 含 子 的 数目 和 分 布 位 置 都 相同 。 内 
含 子 的 变化 也 能 一 定 程 度 地 反映 其 进化 关系 
( Feyereisen, 1999) ,本 研究 为 以 后 对 P450 超 基 
家 族 的 进化 研究 提供 了 基础 。 

CYP6 是 昆虫 特有 的 基因 家 族 , 同 时 也 是 P450 
超 基因 家 族 中 最 主要 的 解毒 酶 基因 家 族 ,对 拟 除 虫 
菊 酯 及 有 机 氧 杀 虫 剂 起 重要 作用 。 不 吉 按 蚊 中 有 3 
个 与 杀 虫 剂 相 关 的 QTL ( quantitative trait locus ) 
(Wondji et al., 2009) , CYP6Y 亚 家 族 基 因 是 其 中 
rp2 QTL 中 的 成 员 (Irving et al., 2012) 4 P] EE TE 
Bic .不 吉 按 蚊 和 微小 按 蚊 的 CYP6P 亚 家 族 中 均 有 基 
能 专 一 代谢 拟 除虫菊 酯 (Miiller et al., 2008; 
Duangkaew et al., 2011; Stevenson et al., 2011; David 
et al., 2013) ,AgCYP6M2 能 代谢 3 类 杀 虫 剂 : 拟 除 忠 
ARS DDT 和 氨基 甲酸 酯 (Mitchell et al., 2012; David 
et al., 2013; Edi et al., 2014)。 因 此 ,本 研究 对 中 华 
fici: CYP6Y ,CYP6M 和 CYP6P 亚 家 族 的 进化 速率 进 
行 分 析 。 相 对 进化 速率 分 析 ( 图 7) 显示 ,中 华 按 蚊 



















































































CYP6Y 和 CYP6M 亚 家 族 基 因 的 进化 速率 显著 快 于 
CYP6P 亚 家 族 , 但 CYP6Y 与 CYP6M 亚 家 族 没有 显著 
差异 ,推测 如 果 中 华 按 蚊 CYP6Y 与 杀 虫 剂 抗 药性 相 
关 , 可 能 与 CYP6M 亚 家 族 相似 可 对 多 类 杀 虫 剂 起 作 
用 。 而 中 华 按 蚊 CYP6Y 是 否 与 杀 虫 剂 抗 性 相关 仍 需 
要 进一步 的 研究 。 

在 冈比亚 按 蚊 吸食 感染 疙 原虫 的 老鼠 血 后 , 差 
异 表达 PASO 可 能 在 细胞 骨架 基因 改变 引起 的 上 皮 
细胞 重 排 中 (在 痉 原 虫 引 起 的 应 激 反 应 或 者 动 合子 
侵 染 和 卵 圳 破裂 过 程 中 ) 起 着 重要 作用 , 其 中 
AgCYP6Y1 和 AgCYP672 在 冈比亚 按 蚊 吸 食 感 染 症 
原虫 的 老鼠 血 1 d 后 ( 即 症 原 虫 动 合子 侵 染 按 蚊 中 
肠 上 皮 细 胞 时 期 ) 分 别 在 脂肪 体 和 中 肠 上 调 表达 
( Félix et al., 2010) 。 而 中 华 按 蚊 中 这 2 个 基因 的 
同 源 基因 AsCYP6Y1 和 AsCYP6Y2 在 症 原虫 感染 中 
华 按 蚊 的 过 程 中 的 作用 以 及 这 2 个 基因 的 功能 还 有 
待 进一步 的 研究 。 


参考 文献 (References) 




















Ai JW, 2008. Genome-Wide Structural and Evolutionary Analysis of 
Cytochrome P450 Genes , and cDNA Cloning and Function of the Key 
Gene CYP18A1 in the Silkworm, Bombyx mori. PhD Dissertation, 
Southwest University, Chongqing. [ WIX, 2008. KÆ (Bombyx 
mori) 全 基因 组 细胞 色素 P450 基因 结构 与 进化 分 析 及 CYPISAT 
克隆 与 功能 研究 . 重庆 : 西南 大 学 博士 学 位 论文 ] 

Ai JW, Wang GH, Li YH, Yu QY, Zhang XS, Zhu Y, Xiang ZH, 


























2008. Molecular cloning, sequence analysis and transcriptional 
activity determination of cytochrome P450 gene CYPI8AI in the 
silkworm, Bombyx mori. Acta Entomologica Sinica, 51(3) : 237 — 
245. [ SEXSISC, 王 根 洪 , FHA, RRA, 张 学 松 , 朱 勇 , 向 仲 
TR, 2008. X P450 基因 CYP18Al 的 克隆 、 序 列 分 析 及 转录 活 
TE. 昆虫 学 报 , 51(3): 237-245] 

Amenya DA, Koekemoer LL, Vaughan A, Morgan J, Brooke BD, Hunt 
RH, Ranson H, Hemingway J, Coetzee M, 2005. Isolation and 























sequence analysis of P450 genes from a pyrethroid resistant colony of 
the major malaria vector Anopheles funestus. Mitochondrial DNA , 16 
(6) : 437 -445. 

David JP, Ismail HM, Chandor-Proust A, Paine MJI, 2013. Role of 
cytochrome P450s in insecticide resistance; impact on the control of 
mosquito-borne diseases and use of insecticides on Earth. 
Philosophical Transactions of the Royal Society B; Biological 
Sciences , 368 (1612) : 20120429. 

Denholm I, Devine GJ, Williamson MS, 2002. Insecticide resistance on 
the move. Science, 297 (5590) : 2222 - 2223. 

Duangkaew P, Pethuan S, Kaewpa D, Boonsuepsakul S, Sarapusit S, 
Rongnoparut P, 2011. Characterization of mosquito CYP6P7 and 
CYP6AA3. differences in substrate preference and kinetic properties. 
Archives of Insect Biochemistry and Physiology, 76(4) : 236 —248. 

Edi CV, Djogbénou L, Jenkins AM, Regna K, Muskavitch MAT, 


672 昆虫 学 报 Acta. Entomologica. Sinica 


57 卷 





Poupardin R, Jones CM, Essandoh J, Kétoh GK, Paine MJI, 
Koudou BG, Donnelly MJ, Ranson H, Weetman D, 2014. CYP6 
P450 enzymes and ACE-1 duplication produce extreme and multiple 
insecticide resistance in the malaria mosquito Anopheles gambiae. 
PLoS Genetics, 10(3) : e1004236. 

Enayati A, Hemingway J, 2010. Malaria management; past, present, 
and future. Annual. Review of Entomology , 55; 569 —591. 

Félix RC, Müller P, Ribeiro V, Ranson H, Silveira H, 2010. 
Plasmodium infection alters Anopheles gambiae detoxification gene 
expression. BMC Genomics , 11(1) : 312. 

Feyereisen R, 1999. Insect P450 enzymes. Annual Review of 
Entomology, 44 (1) : 507 - 533. 

Feyereisen R, 2011. Arthropod CYPomes illustrate the tempo and mode 
in P450 evolution. Biochimica et Biophysica Acta ( BBA) — Proteins 
and Proteomics , 1814 (1) : 19 —28. 

Heckel DG , 2012. Insecticide resistance after silent spring. Science , 337 
(6102) : 1612 —1614. 

Hemingway J, Field L, Vontas J, 2002. An overview of insecticide 
resistance. Science , 298 (5591) ; 96 -97. 

Irving H, Riveron JM, Ibrahim SS, Lobo NF, Wondji CS, 2012. 
Positional cloning of rp2 QTL associates the P450 genes CYP6ZI, 
CYP6Z3 and CYP6M7 with pyrethroid resistance in the malaria 
vector Anopheles funestus. Heredity, 109 (6) ; 383 —392. 

Ji JC, An SH, Li WZ, Luo MH, Yuan GH, Guo XR, 2013. Cloning, 
sequence analysis and prokaryotic expression. of a P450 gene 
HarmCYP9A33 — from 
Noctuidae). Acta Entomologica Sinica, 56(5) ; 465 -474. [Wk 
超 , 安 世 恒 , FAP, 罗 梅 浩 , 原 国 辉 ， 郭 线 站 ,2013. 棉铃 虫 
P450 基因 HarmCYP9433 的 克隆 ,序列 分 析 及 原核 表达 . 昆虫 学 
Jit, 56(5) : 465 - 474] 

Liu QY, Liu XB, 2010. Prevention and control of vector Anopheles: a 


Helicoverpa armigera ( Lepidoptera: 























key approach for malaria elimination in China. Chinese Journal of 
Vector Biology and Control, 21 (5); 409 — 413. [ 刘 起 勇 ， 刘 小 
波 , 2010. 媒介 按 蚊 防 控 : 中 国 症 疾 消除 的 关键 措施 .中 国 媒 
介 生物 学 及 控制 杂志 , 21(5) : 409 - 413] 

Mitchell SN, Stevenson BJ, Miiller P, Wilding CS, Egyir-Yawson A, 
Field SG, Hemingway J, Paine MJI, Randon H, Donnelly MJ, 








2012. Identification and validation of a gene causing cross-resistance 
between insecticide classes in Anopheles gambiae from Ghana. 
Proceedings of the National Academy of Sciences of the United. States 
of America , 109 (16) : 6147 —6152. 

Müller P, Warr E, Stevenson BJ, Pignatelli PM, Morgan JC, Steven A, 
Yawson AE, Mitchell SN, Ranson H, Hemingway J, Paine MJI, 


Donnelly MJ, 2008. Field-caught permethrin-resistant Anopheles 
gambiae overexpress CYP6P3 , a P450 that metabolises pyrethroids. 
PLoS Genetics, 4(11) : e1000286. 

Nelson DR, Koymans L, Kamataki T, Stegeman JJ, Feyereisen R, 
Waxman DJ, Waterman MR, Gotoh O, Coon MJ, Estarbrook RW, 
Gunsalus IC, Nebert DW, 1996. P450 superfamily: update on new 
sequences, gene mapping, accession numbers and nomenclature. 
Pharmacogenetics and Genomics, 6(1) : 1 -42. 

Nikou D, Ranson H, Hemingway J, 2003. An adult-specific CYP6 P450 
gene is overexpressed in a pyrethroid-resistant strain of the malaria 
vector, Anopheles gambiae. Gene, 318; 91 —102. 

Nkya TE, Akhouayri I, Kisinza W, David JP, 2013. Impact of 
environment on mosquito response to pyrethroid insecticides; facts, 
evidences and prospects. /nsect Biochemistry and Molecular Biology , 
43 (4) : 407 - 416. 

Raghavendra K, Niranjan RBP, Prasad GBKS, 2012. Anopheles gambiae 
(Diptera: Culicidae) cytochrome P450 (P450) supergene family; 
phylogenetic analyses and exon-intron organization. Journal of 
Entomology, Ornithology & Herpetology, 1 (1); 1000102. 

Rivero A, Vézilier J, Weill M, Read AF, Gandon S, 2010. Insecticide 
control of vector-borne diseases; when is insecticide resistance a 
problem? PLoS Pathogens, 6(8) : e1001000. 

Stevenson BJ, Bibby J, Pignatelli P, Muangnoicharoen S, O’ Neill PM, 
Lian LY, Müller P, Nikou D, Stenen A, Hemingway J, Sutcliffe 
MJ, Paine MJI, 2011. Cytochrome P450 6M2 from the malaria 
vector Anopheles gambiae metabolizes pyrethroids; sequential 
metabolism of deltamethrin revealed. —/nsect Biochemistry and 
Molecular Biology, 41 (7) : 492 —502. 

Wagner A, 2002. Selection and gene duplication; a view from the 
genome. Genome Biology, 3(5) ; 1012. 

WHO, 2013. World Malaria Report 2013. World Health Organization. 

Wondji CS, Irving H, Morgan J, Lobo NF, Collins FH, Hunt RH, 
Coetzee M, Hemingway J, Ranson H, 2009. Two duplicated P450 
genes are associated with pyrethroid resistance in Anopheles funestus , 
a major malaria vector. Genome Research , 19(3) : 452 —459. 

Yu JF, Hong JF, Wu S, 2012. Resistant status of the four major malaria 
vectors in China. /nternational Journal of Medical Parastic Diseases , 
39(5):314 -317. [ 俞 俊 峰 , WAER, 武松 , 2012. 我 国 四 大 传 

疙 媒介 抗 性 现状 .国际 医学 寄生 虫 病 杂 志 , 39(5) : 314-317] 

Zhou Q, Wang W, 2004. Detecting natural selection at the DNA level. 
Zoological Research, 25(1); 73 —80. [ 周 琦 , XE X, 2004. DNA 
水 平 自然 选择 作用 的 检测 .动物 学 研究 , 25(1) : 73 -80] 

(责任 编辑 : 赵 利 辉 ) 




















